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ABSTRACT 


Free  electron  lasers  (FEL's)  are  devices  used  worldwide 
for  several  purposes.  In  the  military,  especially  in  the 
Navy,  they  can  be  used  for  self-defense  against  missiles, 
and  small  boats.  Installed  on  a  ship,  an  FEL  represents  a 
multi-mission,  deep  magazine,  long  range  weapon.  This  thesis 
will  describe  briefly  the  basic  components  and  principles  of 
operation.  It  also  explores,  by  simulations,  the  effects  of 
changing  some  of  the  parameters  that  generate  the  laser 
beam. 


V 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


VI 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . 1 

A.  HISTORY  OF  THE  LASER . 1 

B.  HISTORY  OF  A  FREE  ELECTRON  LASER . 2 

C.  ADVANTAGES  AND  DISADVANTAGES  OF  THE  FREE  ELECTRON 

LASER . 3 

II.  BASIC  COMPONENTS  OF  A  FREE  ELECTRON  LASER . 5 

A.  THE  INJECTOR . 6 

B .  LINEAR  ACCELERATOR . 7 

C .  UNDULATOR . 8 

D .  RESONATOR . 9 

E .  SEED  LASER . 9 

F.  ACCESSORIES . 10 

G.  ELECTRON  BEAM  DUMP . 10 

III.  BASIC  THEORY  OF  A  FREE  ELECTRON  LASER . 11 

A.  THE  PENDULUM  EQUATION . 11 

B.  PHASE  SPACE . 13 

C.  ELECTRON  BEAM  AND  OPTICAL  BEAM  INTERACTION . 15 

D.  FEL  WAVE  EQUATION . 16 

E.  WEAK  FIELD  GAIN . 17 

1.  Low  Gain  Regime  (j  ^  n) . 18 

2.  High  Gain  Regime  (j»n) . 19 

F.  STRONG  FIELD  GAIN . 19 

G.  TAPERED  WIGGLER . 20 

1.  Linear  Tapering . 21 

2 .  Step-Tapering . 22 

H.  DIFFRACTION . 26 

I.  SHORT  PULSES  AND  SLIPPAGE . 28 

IV.  SIMULATIONS  OF  FREE  ELECTRON  LASERS  IN  AMPLIFIER 

CONFIGURATION . 31 

A.  DESCRIPTION  OF  THE  INITIAL  PARAMETERS . 31 

B.  SIMULATIONS  OF  THE  FIRST  FEL . 33 

C.  SIMULATIONS  OF  THE  SECOND  FEL . 46 

D.  SIMULATIONS  OF  THE  THIRD  FEL . 54 

V.  CONCLUSION . 61 

LIST  OF  REFERENCES . 63 

INITIAL  DISTRIBUTION  LIST . 65 


VI 1 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


VI 11 


LIST  OF  FIGURES 


Figure  1.  Basic  components  of  a  free  electron  laser 

(oscillator  mode)  from  [3] . 5 

Figure  2.  Basic  scheme  of  a  photocathode  injector  from 

[4] . 6 

Figure  3.  Schematic  of  a  linear  accelerator  from  [4] . 7 

Figure  4.  Schematic  of  an  undulator  from  [6] . 8 

Figure  5.  Schematic  of  a  resonator  from  [3] . 9 

Figure  6.  Arrangement  of  a  seed  laser  from  [5] . 10 

Figure  7 .  Separatrix  and  orbits  on  a  pendulum  phase  space 

diagram . 14 

Figure  8.  FEL  phase  space  diagram  from  [3] . 15 

Figure  9.  FEL  electron's  bunching  from  [3] . 16 

Eigure  10.  Positive  linear  tapered  undulator  from  [3] . 21 

Eigure  11.  Negative  linear  tapered  undulator  from  [3] . 22 

Eigure  12.  Positive  step-tapered  undulator  from  [3] . 22 

Eigure  13.  Negative  step-tapered  undulator  from  [3] . 22 

Eigure  14.  Phase  space  diagram  of  a  positive  linearly 

tapered  undulator  from  [3] . 24 

Eigure  15.  Phase  space  diagram  of  a  negative  linearly 

tapered  undulator  from  [3] . 25 

Eigure  16.  Phase  space  diagram  of  a  positive  step-tapered 

undulator  from  [3] . 25 

Eigure  17.  Phase  space  diagram  of  a  negative  step-tapered 

undulator  from  [3] . 26 

Eigure  18.  Electron  and  light  slippage  from  [3] . 29 

Eigure  19.  Optical  lethargy  from  [3] . 30 

Eigure  20.  Graph  of  extraction  vs.  taper  start  times . 35 

Eigure  21.  Eull  simulation  output  for  taper  start  time 

0.7 . 36 

Eigure  22.  Graph  of  extraction  vs.  linear  taper  rates . 37 

Eigure  23.  Graph  of  extraction  vs.  linear  taper  rates  for 

N=14  0 . 3  9 

Eigure  24.  Electron  bunching  and  optical  field  for  linear 

taper  rate  5=30n . 39 

Eigure  25.  Graph  of  extraction  vs.  electron  beam  pulse 

duration . 40 

Eigure  26.  Extraction  vs.  dimensionless  taper  start  time 

for  micropulse  charge  of  0.4  nC . 42 

Eigure  27.  Extraction  vs.  dimensionless  taper  rates  for 

micropulse  charge  of  0.4  nC . 43 

Eigure  28.  Results  for  micropulse  charge  of  0.4  nC  and 

taper  rate  5=60n . 44 


IX 


Figure  29.  Extraction  vs.  dimensionless  taper  rates  for 

emittance=5 . 1mm  mrad . 46 

Figure  30.  Extraction  vs.  number  of  periods  of  the  second 

EEL . 4  9 

Figure  31.  Results  for  N=140  of  the  second  EEL . 50 

Figure  32.  Extraction  vs.  dimensionless  taper  start  time 

of  the  second  FEL . 51 

Figure  33.  Extraction  vs.  dimensionless  linear  taper  rates 

of  the  second  FEL . 52 

Figure  34.  Results  for  5  =  lOOn  of  the  second  FEL . 53 

Figure  35.  Extraction  vs.  number  of  periods  of  the  third 

FEL . 5  6 

Figure  36.  Extraction  vs.  dimensionless  taper  start  time 

of  the  third  FEL . 57 

Figure  37.  Extraction  vs.  dimensionless  taper  rates  of  the 

third  FEL . 58 

Figure  38.  Results  for  N=130  of  the  third  FEL . 59 


X 


LIST  OF  TABLES 


Table  1 .  Undulator  parameters  used  to  describe  an  FEL  in 

a  simulation . 32 

Table  2 .  Electron  beam  parameters  used  to  describe  an 

EEL  in  a  simulation . 32 

Table  3.  Optical  parameters  used  to  describe  an  FEL  in  a 

simulation . 33 

Table  4.  Initial  electron  beam  parameters  of  the  first 

simulated  FEL  device . 33 

Table  5.  Initial  undulator  parameters  of  the  first 

simulated  EEL  device . 34 

Table  6.  Initial  optical  parameters  of  the  first 

simulated  EEL  device . 34 

Table  7.  Dimensionless  parameters  of  the  first  simulated 

FEL  device . 34 

Table  8.  New  parameters  for  N=140 . 38 

Table  9.  Parameters  that  changed  for  micropulse  charge 

of  0.4 . 41 

Table  10.  Parameters  that  changed  after  changing  the 

emittance . 45 

Table  11.  Initial  electron  beam  parameters  of  the  second 

simulated  FEL . 47 

Table  12.  Initial  undulator  parameters  of  the  second 

simulated  EEL  device . 47 

Table  13.  Initial  optical  parameters  of  the  second 

simulated  EEL  device . 48 

Table  14.  Dimensionless  parameters  of  the  second 

simulated  EEL  device . 48 

Table  15.  Initial  electron  beam  parameters  of  the  third 

simulated  EEL . 54 

Table  16.  Initial  undulator  parameters  of  the  third 

simulated  EEL  device . 55 

Table  17.  Initial  optical  parameters  of  the  third 

simulated  FEL  device . 55 

Table  18.  Dimensionless  parameters  of  the  third  simulated 

FEL  device . 55 


XI 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


Xll 


ACKNOWLEDGMENTS 


I'd  like  to  thank  Professors  Bill  Colson,  Joe  Blau  and 
Robert  Armstead  for  their  aid  in  the  preparation  of  this 
work,  without  them  this  work  would  not  be  possible. 

Special  thanks  to  my  wife,  Rosa,  for  her  understanding 
and  help  during  the  time  I  spent  in  this  work. 


Xlll 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


XIV 


I. 


INTRODUCTION 


A.  HISTORY  OF  THE  LASER 

The  word  "laser"  is  an  acronym  for  "light  amplification 
by  the  stimulated  emission  of  radiation".  Laser  is  an 
optical  source  that  emits  nearly  monochromatic  light  in  a 
narrow  coherent  beam. 

The  laser  was  proposed  as  a  variation  of  the  maser 
("Microwave  amplification  by  the  stimulated  emission  of 
radiation")  principle  in  the  late  1950s  and  was  demonstrated 
in  July  1960  by  Theodore  Maiman  at  Hughes  research 
laboratories  [1]. 

Albert  Einstein,  in  his  1916  paper,  Strahlungs-Emission 
und-ahsorption  nach  der  Quantentheorie ,  laid  the  foundation 
of  the  laser  and  the  maser.  He  applied  Max  Planck's  law  of 
radiation  in  terms  of  the  probability  of  absorption  and 
spontaneous  and  stimulated  emission.  In  1928,  Rudolph  W. 
Landenburg  confirmed  the  existence  of  stimulated  emission 
and  negative  absorption  experimentally.  In  1939,  the 
Russian,  Valentin  A.  Fabrikant,  suggested  the  use  of 
stimulated  emission  to  amplify  "Short"  waves.  In  1947, 
Willis  E.  Lamb  and  R.C.  Retherford  demonstrated  apparent 
stimulated  emission  in  hydrogen  spectra.  Finally  in  1950 
Alfred  Kastler  proposed  the  method  of  optical  pumping,  which 
was  experimentally  confirmed  by  Brossel,  Kastler  and  Winter 
two  years  later  [1] . 

In  1953,  Charles  H  Townes,  James  P  Gordon  and  Herbert 
J.  Zeiger  made  the  first  "Maser"  to  produce  microwave 
radiation.  In  1957,  Townes  and  Arthur  L.  Schawlow  began  to 
produce  emissions  in  the  infrared  band  and  then  adapted  it 
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to  optical  wavelengths.  In  1958,  Prokhorov  proposed  the 
usage  of  an  open  resonator  in  the  maser,  and  Schawlow  and 
Townes  settled  on  an  open  resonator  design [1] . 

The  first  working  laser  was  made  by  Theodore  Maiman  in 
1960  in  Malibu,  California;  a  collaboration  of  several 
institutions.  This  laser  was  made  with  a  synthetic  ruby 
crystal,  pumped  by  a  flashlamp  to  produce  light  at 
wavelength  694  nm  [1] . 

In  1960,  All  Javan,  William  Bennet  and  Donald  Heriot 
made  the  first  gas  laser  using  Helium  and  Neon  [1] . 

The  first  laser  diode  was  made  by  Robert  N.  Hall  in 
1962.  It  was  made  of  gallium  arsenide  and  emitted  in  the 
near-infrared  at  850  nm.  The  first  semiconductor  laser  with 
visible  emission  was  made  the  same  year  by  Nick  Holonyak, 
Jr.  These  semiconductor  lasers  could  be  used  only  in  pulsed 
operation,  and  indeed  only  when  cooled  to  liquid  nitrogen 
temperatures  (77  K)  [1] . 

In  1970,  Zhores  Alferov  in  the  Soviet  Union  and  Izuo 
Hayashi  and  Morton  Panish  of  Bell  Telephone  Laboratories 
independently  developed  laser  diodes  continuously  operating 
at  room  temperature,  using  the  heteroj unction  structure  [1] . 

B.  HISTORY  OF  A  FREE  ELECTRON  LASER 

The  Free  Electron  Laser  (EEL)  appeared  in  the  middle 
1970s,  with  the  well-documented  experiments  at  Stanford  by 
Dr.  John  M.J.  Madey  (professor  at  that  university  until 
1988)  .  The  EEL  by  itself  has  a  background  which  dates  back 
to  at  least  1947  and  a  patent  by  Elmer  Corn,  in  which  he 
describes  amplification  producing  fast  electromagnetic  wave 
(periodic  electron  beam)  interactions  [2]. 
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In  the  1950s,  Motz  worked  on  an  undulator,  and  he 
passed  a  beam  from  a  3-MeV  accelerator  through  it  to  produce 
millimeter  wave  radiation.  He  later  used  a  100-MeV  beam  to 
produce  visible  light.  Another  experiment  of  Motz  was  the 
Ubitron,  a  mildly  relativistic  microwave  FEL.  The  power 
capability  of  this  device  was  so  much  greater  than  that  of 
conventional  klystrons  and  magnetrons  that  records  for  peak 
power  generation,  which  would  stand  for  two  decades,  were 
established  at  both  centimeter  and  millimeter  wavelengths 
[2]  . 

Nowadays  the  tendency  of  the  institutions  is  to  design 
EELS  to  operate  at  shorter  wavelengths.  The  most  powerful 
EEL  in  the  United  States  is  located  in  Jefferson' s  Lab  at 
Newport  News,  Virginia;  having  a  power  of  14  kilowatts,  this 
device  works  in  an  oscillator  configuration. 

C.  ADVANTAGES  AND  DISADVANTAGES  OF  THE  FREE  ELECTRON 

LASER 

The  Free  Electron  laser  has  several  advantages  over 
other  lasers.  It  is  more  flexible  than  the  others;  for 
example,  it  is  tunable  by  design.  That  means  that  the  output 
wavelength  can  be  changed  as  desired  by  changing  the  energy 
of  the  electron  beam  or  the  properties  of  the  undulator.  A 
second  advantage  is  that  the  EEL  can  operate  continuously 
(24  hours  a  day)  at  a  relatively  low  cost,  whereas  chemical 
lasers  have  to  be  replenished  after  a  short  time  of 
operation,  and  they  may  generate  toxic  exhausts.  The  EEL 
will  operate  all  the  time  if  it  has  a  continuous  power 
supply  and  will  not  damage  the  environment  with  a  toxic 
exhaust.  Also,  EEL  can  reach  higher  powers  than  the  chemical 
lasers,  because  their  lasing  medium  does  not  get  hot. 
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II.  BASIC  COMPONENTS  OF  A  FREE  ELECTRON  LASER 


There  are  two  basic  configurations  of  a  free  electron 
laser:  The  oscillator  configuration  and  the  amplifier 
configuration.  Both  configurations  have  many  common 
components;  the  difference  is  in  the  way  that  the  laser  beam 
is  generated.  The  basic  components  are:  injector,  linear 
accelerator,  undulator,  resonator  (or,  for  the  amplifier,  a 
laser  seed),  beam  dump  and  accessories.  Figure  1  shows  the 
basic  components  of  a  free  electron  laser. 


Electron  Accelerator 


Figure  1.  Basic  components  of  a  free  electron  laser 

(oscillator  mode)  from  [3]. 
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A. 


THE  INJECTOR 


This  is  the  first  part  of  the  device  and  is  the  source 
of  the  electron  beam.  Basically  it  contains  a  cathode,  as 
does  a  TV  tube.  The  electron  beam  is  generated  in  short 
pulses  (on  the  order  of  picoseconds)  with  energy  of  5  MeV 
approximately,  a  value  common  to  both  configurations.  There 
are  several  types  of  injectors;  among  these  we  can  find  the 
direct-current  injector  (DC),  the  radio-frequency  injector 
(RF)  and  the  semiconducting  radio-frequency  injector;  and 
all  of  them  can  be  used  in  an  FEL  system.  Some  injectors  are 
driven  by  a  laser  in  order  to  provide  the  required  beam;  the 
drive  laser  consists  of  a  diode-pumped  NdiGlass  oscillator 
and  a  chirped  pulse  amplification  system  consisting  of  a 
granting  stretcher,  a  f lashlamp-pumped  Nd: Glass  regenerative 
amplifier  and  a  grating  compressor  [8]  .  Figure  2  shows  a 
basic  scheme  of  a  photocathode  injector. 


Figure  2 . 


Basic  scheme  of  a  photocathode  injector  from 

[4]  . 
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B. 


LINEAR  ACCELERATOR 


The  linear  accelerator  (linac)  is  a  critical  component 
of  the  equipment,  and  increases  the  energy  of  the  electron 
beam  from  5  MeV  at  the  output  of  the  injector  to  about  100 
MeV.  This  increment  of  energy  is  done  in  the  accelerator  RF 
cavities.  Electrons  gain  energy  from  the  RF  fields  as  they 
go  from  one  cavity  to  the  next,  traveling  through  an 
electric  field.  The  efficiency  and  gradient  of  the  linac, 
and  therefore  the  entire  system,  can  be  increased  by  use  of 
a  superconducting  accelerator,  which  must  operate  at 
cryogenic  temperatures.  This  requires  the  installation  of  a 
liquid  helium  refrigerator.  Figure  3  shows  the  schematic  of 
a  linear  accelerator. 


Seven  cell  cavity  is  incorporated  in  a  devvar. 


cavities  make  a  cryomodule 


Figure  3.  Schematic  of  a  linear  accelerator  from  [4] . 
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c. 


UNDULATOR 


The  undulator  is  where  the  laser  beam  is  created  and 
enhanced,  and  where  much  of  the  interesting  physics  occur. 
An  undulator  is  a  series  of  alternating  magnets,  typically 
forming  about  20  -  100  periods,  separated  by  a  few 

centimeters,  and  delivering  a  typical  magnetic  peak  field  of 
about  1  Tesla.  The  characteristics  of  the  laser  beam  are 
determined  by  the  undulator  design,  being  a  region  in  which 
the  electrons  are  made  to  wiggle  by  a  static  alternating 
magnetic  field.  Since  wiggling  electrons  radiate  light,  we 
can  store  this  light  in  a  resonator;  and  we  can  stimulate 
emission  by  radiating  in  the  presence  of  this  same  light. 
While  electrons  are  in  the  presence  of  the  light  there  is  an 
exchange  of  energy  to  and  from  the  light.  The  electron 
distribution  in  phase  space  must  be  arranged  in  order  to 
obtain  the  most  efficient  transfer  of  energy  from  the 
electrons  to  the  light.  Additional  information  will  be 
written  about  this  effect  in  a  later  chapter.  In  the 
amplifier  configuration,  the  light  passes  only  one  time 
through  the  undulator,  and  then  we  require  a  longer 
undulator  to  extract  the  largest  possible  energy  from  the 
electron  beam.  Figure  4  shows  the  basic  architecture  of  an 
undulator . 


Figure  4.  Schematic  of  an  undulator  from  [6] . 
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D. 


RESONATOR 


The  resonator  is  the  component  where  the  light  is 
stored  in  the  oscillator  configuration.  It  consists  of  an 
evacuated  cavity  bracketed  by  two  mirrors  typically 

separated  by  20  -  30  meters,  between  which  the  light  is 

stored.  One  of  the  mirrors  is  perfectly  reflective  (100%), 

the  other  is  partially  transmitting  (it  transmits  some 
percentage  of  the  light) .  It  is  through  this  mirror  that  the 
laser  beam  escapes  to  the  outside.  Resonators  are  often 
characterized  by  their  quality  factor  Q,  which  is 
established  by  the  inverse  of  their  loss  per  pass  (1  /  Loss 
per  pass) .  Figure  5  shows  a  basic  arrangement  of  a 

resonator . 


electron  beam  pulses 
^  wiggler 


100  % 

reflecting 

mirror 


•  •> 

•  ^ 

partially 

transmitting 

mirror 


optical 

laser 

pulses 


Figure  5  . 


Schematic 


of  a  resonator  from  [3]  . 


E .  SEED  LASER 

For  an  amplifier  configuration  it  is  necessary  to  have 
an  external  laser  source,  because  unlike  the  oscillator,  no 
light  is  stored  in  the  optical  cavity.  This  external  laser 
source  is  the  "seed  laser";  it  generates  a  laser  beam  of 
weak  field  (low  power  usually  10  -  100  W)  that  is  introduced 
to  the  undulator  for  amplification  to  the  desired  output  (1 
MW  in  our  case)  .  Figure  6  shows  a  basic  arrangement  of  a 
seed  laser. 
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High-Power 
PEL  Beam 


Figure  6.  Arrangement  of  a  seed  laser  from  [5] . 

F .  ACCESSORIES 

"Accessories"  implies  all  the  devices  that  control  the 
behavior  of  the  beams  (electron  and  optical) .  They  consist 
of  bending  and  focusing  magnets  that  center  the  electron 
beam  and  change  its  direction,  and  the  mirrors  that  direct 
the  optical  beam. 

G.  ELECTRON  BEAM  DUMP 

It  is  the  last  part  of  the  electron  path.  In  an  FEL 
with  recirculation  of  the  electron  beam,  the  electrons  enter 
to  the  linac  a  second  time  180°  out  of  phase  with  the  RF 
fields,  giving  most  of  their  energy  back  to  the  electric 
field.  At  the  output  of  the  linac,  the  recirculated  electron 
beam  has  energy  of  about  5  MeV,  which  is  about  the  same 
energy  that  they  have  when  they  are  generated  at  the 
injector.  These  electrons,  with  the  residual  energy,  are 
dumped  into  what  is  called  the  "electron  beam  dump."  At  that 
relatively  low  energy,  the  electrons  create  less  heat  and 
radiation,  so  less  shielding  is  required  on  this  electron 
beam  dump . 
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III.  BASIC  THEORY  OF  A  FREE  ELECTRON  LASER 


A.  THE  PENDULUM  EQUATION 


Since  an  accelerated  charged  particle  emits  radiation, 
we  can  create  radiation  by  accelerating  electrons  inside  the 
undulator.  We  do  this  by  making  the  relativistic  electrons 
'wiggle'  through  the  undulator.  When  an  electron  enters  to 
the  undulator,  it  goes  through  electric  and  magnetic  fields 
and  encounters  a  force  that  acts  upon  it  given  by  the 
Lorentz  force  equation.  For  relativistic  electrons  in  the 
undulator,  the  complete  Lorentz  FEL  force  equations  are 
given  by  [3] 


d[rp) 

dt 


dy 

dt 


-  —  pXE, 
me 


(III.l) 

(III  .2) 


\  =  l-p\  (III. 3) 

r 

where  e  is  the  electron  charge  magnitude,  m  the  electron 
mass,  E  the  electric  field,  B  the  magnetic  field,  y0  =  v/c. 


V  is  the  electron  velocity,  c  is  the  speed  of  light  and  y 
is  the  relativistic  Lorentz  factor. 

If  we  substitute  the  electric  and  magnetic  fields  of 
the  FEL  into  equation  (III.l),  we  get 


d(r/^±) 

dt 


d[rP,) 

dt 


e 

me 


E[\-  P^)[cosy/,-smy/,Qi)  + 


e 

me 


E[l3,co^y/  -  +  B[/3^?,m{k^z)- 


,cos(^„z),0)] 

cos(^„z)) 


(III. 4) 


(III . 5) 
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where  +lj  is  the  magnitude  of  the  longitudinal 

z  component  of  the  electron  velocity,  E  is  the  magnitude  of 
the  electric  field,  B  is  the  magnitude  of  the  magnetic  field 
and  xj/  is  the  phase.  We  can  solve  for  the  transverse 

velocity  components  = -(ii/7)(cos(A:^z),sin(A:„z),0)  where  fi^  and 

(5^,  are  the  components  of  the  velocity  and  K  =  eBX^Ilnmc^  . 

Substituting  the  transverse  velocity  components 
into  (III. 2)  we  get 

°  e  -  -  cKE  cKE 

y  = - P'E  = - (cosA:_,zcos^y-sinA:^zsin^)  = - cos (4"  +  ^) ,  ( III .  6) 

me  ymc  ymc 

where  P  =  z  -  cot  is  the  electron  phase,  which  represents 

the  position  of  the  electron  within  an  optical  wavelength. 
If  we  take  the  derivative  of  the  electron  phase  with  respect 
to  the  dimensionless  time  z^etjE,  where  L  =  NX^  is  the  length 

of  the  undulator,  and  N  the  number  of  periods  of  the 
undulator,  we  get  the  electron  phase  velocity 


v  =  l  =  L\{k  +  kP)p^-k']  . 

For  highly  relativistic  electrons  (f»l),  we  can  approximate 
,  and  near  resonance  [k»kP^,  so  that  the 
derivative  of  the  electron  phase  becomes 


P  =  L  [k  +  kP)(\-[K^  +\)l2y^Y^  =Lk(\-[K^ +\)l2y^Y 


If  we  take  the  second  derivative  of  the  electron  phase,  we 
get 


yPk 
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Substituting  (III.  6)  into  the  last  equation  and  defining 
\a\  =  A7T:NeK\E\L j Y^mc^  as  the  dimensionless  field  amplitude,  we 
obtain  the  dimensionless  pendulum  equation, 

=  K  =  |a|cos(i^  +  ^z))  .  (III. 7) 

This  equation  describes  the  microscopic  motion  of  electrons 
in  phase  space  within  the  undulator. 

B .  PHASE  SPACE 

The  phase  space  representation  is  a  diagram  that  shows 
the  evolution  in  position  and  velocity  coordinates.  It  is 
useful  to  show,  as  in  the  last  section,  the  electron  energy 
exchange  in  an  FEL.  This  diagram  shows  the  electron  phase 
velocity  versus  the  phase.  There  are  two  types  of  orbits  on 
a  phase  space  diagram:  closed  and  open  orbits. 

Closed  orbits  for  a  classical  pendulum  represent  a 
pendulum  swinging  back  and  forth;  open  orbits  occur  when  the 
pendulum  goes  over  the  top  due  to  a  large  angular  velocity. 
The  border  between  those  two  states  is  called  the 
separatrix.  Figure  7  shows  a  separatrix  (black  line)  and  the 
closed  (green  line)  and  open  (red  doted  line)  orbits  on  a 
phase  space  diagram;  horizontal  axis  is  the  angular 
displacement  and  vertical  axis  the  angular  velocity. 
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Figure  7 .  Separatrix  and  orbits  on  a  pendulum  phase 

space  diagram. 

Since  we  suggested  the  similarity  between  the  FEL  and  a 
pendulum,  and  we  obtained  a  pendulum  equation  describing  the 
behavior  of  electrons  in  an  FEL,  we  can  state  that  |a|  and  (j) 

are  roughly  constant  for  the  case  of  low  gain  in  the  FEL;  in 
the  same  way,  the  electron  phase  4"  can  be  viewed  as  the 
electron' s  position  within  an  optical  wavelength,  and  v  is 
its  phase  velocity  relative  to  resonance  (v=0) . 

For  the  FEL,  the  separatrix  is  given  by  [3] : 

=  2|a|[l  +  sin(4'  +  ^zJ)]  (III. 8) 

Figure  8  shows  a  phase  space  diagram  of  the  behavior  of 
approximately  20  electrons  evenly  spaced  in  ^  with  vo=0  in  a 
free  electron  laser.  Yellow  color  indicates  the  electrons  at 
the  beginning  of  the  undulator  (t=0)  and  blue  color 
indicates  the  electrons  at  the  end  of  the  undulator  (t=1)  . 
In  this  particular  situation  the  gain  of  the  system  is  zero, 
because  half  of  the  electrons  gain  energy  and  half  of  the 
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electrons  lose  it.  In  order  to  get  gain,  it  is  necessary  to 
initiate  the  process  with  the  electrons  off  resonance  as 
Figure  9  shows . 


5 


V 


-5 


Figure  8.  FEL  phase  space  diagram  from  [3] 

C.  ELECTRON  BEAM  AND  OPTICAL  BEAM  INTERACTION 

One  important  issue  in  an  FEL  is  the  energy  exchange 
between  electrons  and  the  optical  field.  As  noted 
previously,  values  of  v  close  to  resonance  (v=0)  are 
required  for  optimal  energy  exchange.  Phase  space  evolutions 
of  the  electrons  showing  a  good  bunching  or  grouping  result 
in  bigger  gain.  So  it  is  necessary  to  create  an  optimal 
bunching  to  get  the  best  values  of  gain. 

Figure  9  shows  a  phase  space  diagram  of  the  evolution 
of  the  electrons  in  a  free  electron  laser  for  a  more 
realistic  beam  with  1000  particles  and  a  small  initial 
energy  spread. 
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Figure  9.  FEL  electron's  bunching  from  [3] . 


D.  FEL  WAVE  EQUATION 


From  Maxwell's  equations,  we  can  obtain  a  full  wave 
equation  for  the  FEL  optical  field  in  terms  of  the 
transverse  current  [3] : 


2 ,  1  5  A  An 

V  A — = - , 

c  ot  c 


(III . 9) 


where  A  is  the  optical  magnetic  vector  potential  and  J  is 
the  current  density. 

For  our  helical  undulator,  a  form  for  the  optical 
vector  potential  can  be  written  as  [3] : 


(III .  10) 


where  E(z,t)  is  the  complex  optical  amplitude,  a  =  kz-cot  is 
the  carrier  wave  phase,  i  =  (-/,l,0)  is  the  optical  field' s 
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polarization  vector  for  circularly  polarized  light,  and 
X  =  xi +  _yj  +  zk  is  the  position  vector. 

Substituting  the  solution  in  the  wave  equation  we  can 
get  the  results  for  the  four  dimensions.  We  can  make  some 
simplifications  and  substitutions  like  the  slowly  varying 
amplitude  and  phase  approximation  that  will  be  discussed 
later  in  this  chapter,  in  order  to  get  the  FEL  wave  equation 
as  [  3  ]  : 

T 

dr 

where  a(r)  =  |a|e'^  is  the  complex  dimensionless  optical  field 
amplitude;  the  dimensionless  current  is  j  =  'iN{e7T:KLY  n^j y^mc^  , 
where  is  the  beam  particle  density,  measuring  the 

interaction  between  the  electron  beam  and  the  optical  mode; 
and  <e^“"  >  measures  the  amount  of  electron  bunching.  If  j  is 
large  (j  >>  n)  ,  the  optical  field  a  changes  rapidly,  so  we 
have  high  gain.  If  j  is  small  (j  ^  n) ,  we  have  low  gain. 


)_ 


<  -ye > 


(III . 11) 


E. 


WEAK  FIELD  GAIN 


As  mentioned  previously,  we  can  obtain  gain  through  the 
interaction  of  the  bunched  electron  beam  with  the  optical 
field  inside  the  undulator.  An  optical  field  is  considered 
weak  when  |a|<;r  it  is  smaller  than  n.  The  following  wave 
equations  are  valid  for  weak  and  strong  optical  fields  [3] : 

dial 


dr 


=  -  j  (cosi^^  +  (f)'^  ,  and 


(III . 12) 


^  =  A(sin(^  +  ^))  . 


(Ill .  13) 
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They  describe  the  change  in  optical  amplitude  |a|  and 

phase  cp  during  FEL  interaction  with  the  electron  beam.  The 
pointing  brackets  denote  the  average  of  all  the  electrons 
sampled  within  an  optical  wavelength  of  the  beam. 

From  the  coupling  between  the  electron  beam  and  optical 
field,  we  can  calculate  the  dimensionless  current  density  j. 


1.  Low  Gain  Regime  (j  ^  n) 


In  this  regime,  since  j  is  small  there  is  not  an 
appreciable  change  in  the  optical  field  amplitude  or  phase 
as  the  electrons  pass  through  the  undulator.  The  evolution 
of  the  optical  field  amplitude  and  phase  versus  i 
(dimensionless  time  for  the  electrons  to  pass  through  the 
undulator  from  t=0  through  t=1)  can  be  obtained  by  expanding 
^  and  V  in  powers  of  |a|  =  ao  to  get  [7] 


!(r)  =, 


1  +  7 


2-2cos(Kor)-KQrsin(Kor) 

2^1 


(III . 14) 


2sin(Kor)-KQr(l  +  cos(KQr)) 


2k' 


(III . 15) 


where  ag  is  the  initial  optical  field  amplitude,  and  Kq  is 
the  initial  phase  velocity  of  the  electron  beam. 

In  an  FEL  system,  the  gain  is  sensitive  to  the  initial 
phase  velocity  Vq  .  If  the  electron  beam  starts  with  Kq  =  0 
(at  resonance),  the  resultant  gain  is  negligible.  The 
resulting  gain  inside  the  undulator  assuming  weak  fields  is 
given  by  [7] 


G{^)  =  J 


2  -  2  cos  ( Kgr )  -  Kgr  sin  ( Kgr ) 


(III . 16) 
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2.  High  Gain  Regime  (j»n) 

In  this  regime,  the  dimensionless  current  density  is 
large,  and  it  strongly  affects  the  optical  field  amplitude 
and  phase  as  the  electrons  go  through  the  undulator.  In  this 
case,  the  optical  field  amplitude  and  phase  during  the 
electron  beam  interaction  with  the  optical  field  are  given 
by  [7] 

|a(r)|  =  — ^  (III. 17) 

I  V  ;i  3 


(III .  18) 


From  these  formulas,  we  can  see  large  changes  in  the 
amplitude  and  phase  due  to  the  current  density;  then  for 
large  changes  in  the  phase,  the  separatrix  is  shifted  and 
the  optimum  bunching  location  changes  on  the  phase  space 
diagram.  In  contrast  to  the  low  gain  regime,  the  gain  is  not 
as  sensitive  to  the  initial  phase  velocity,  but  the  gain 
spectrum  changes  as  well. 

The  expression  for  the  gain  in  the  undulator  for  this 
regime  is  given  by  [7] 


(III .  19) 


F.  STRONG  FIELD  GAIN 

For  this  case,  in  both  regimes  (low  and  high  gain)  , 
strong  optical  fields  make  the  electron  beam  overbunch  in 
the  undulator  as  the  electrons  travel  in  it.  This  process  is 
called  saturation.  Overbunching  creates  absorption  of  energy 
by  the  electrons  from  the  optical  field,  reducing  the  gain 
of  the  system.  In  the  low  gain  regime,  saturation  is 
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achieved  at  much  lower  initial  field  amplitudes  than 

the  high  gain  regime  « 2(7/2)^^^  [7]. 

When  there  is  not  saturation,  the  phase  space  evolution 
of  the  electrons  is  slower  and  the  bunch  of  electrons  travel 
downward  to  a  phase  of  ^  ~  n.  However,  the  stronger  optical 
field  causes  a  greater  force  on  the  electrons  and  more 
evolution  in  phase  space.  The  bunch  forms  and  continues 
along  closed  orbits  until  it  begins  an  upward  trend,  taking 
the  energy  away  from  the  optical  field. 

G.  TAPERED  WIGGLER 

In  this  work,  an  analysis  of  an  FEL  in  the  amplifier 
configuration  is  made;  this  type  of  system  requires  a  much 
longer  undulator  than  the  oscillator.  If  we  want  to  increase 
the  gain  and  extraction  of  the  system,  we  can  extend  the 
saturation  of  electrons  by  tapering  the  undulator.  Taper  is 

Pi 

a  change  in  the  resonance  condition  X  =  +  as  the 

electrons  travel  along  the  undulator;  this  alteration 
extends  the  energy  exchange  between  electrons  and  the 
optical  field;  this  extended  exchange  allows  a  greater 
amount  of  energy  taken  from  the  electrons  and  an  increased 
extraction  of  energy  from  the  device.  As  the  electrons 
evolve,  they  loose  energy,  and  y  decreases,  then,  in  order 
to  keep  the  same  wavelength,  or  K  must  be  decreased. 

There  are  several  ways  to  taper  the  undulator;  the  most 
common  is  to  vary  the  undulator  parameter  K  by  changing  the 
magnetic  fields  along  it.  The  easiest  way  to  do  this  taper 
is  to  change  the  size  of  the  gap  between  magnets  by  a  small 
amount  (Tenths  of  millimeters)  ;  this  will  change  the 
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strength  of  the  field  at  certain  points.  In  this  way,  we  can 
generate  different  types  of  tapering. 

1 .  Linear  Tapering 

This  tapering  method  consists  of  a  continuous  increment 
or  decrement  of  the  magnetic  field  strength  along  the 
undulator.  This  variation  of  the  field  can  be  negative 
(increasing  gap)  or  positive  (decreasing  gap) .  The  taper  can 
start  at  any  location  along  the  undulator  and  continues 
linearly  through  the  rest  of  the  undulator.  Typically,  the 
taper  stars  near  the  location  where  the  optical  field  is 
expected  to  reach  its  saturation  value  . 

As  we  can  see  in  Figures  10  and  11  (positive  and 
negative  linear  taper) ,  this  technique  exposes  the  electron 
beam  to  a  linearly  changing  magnetic  field. 


Figure  10. 


Positive  linear  tapered  undulator  from  [3] . 
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Figure  11.  Negative  linear  tapered  undulator  from  [3] . 

2 .  Step-Tapering 

Other  useful  taper  method  is  the  step  taper.  This 
method  consists  of  changing  the  gap  at  a  certain  point 
within  the  undulator  and  maintaining  this  gap  through  the 
rest  of  it.  Figures  12  and  13  show  a  positive  and  a  negative 
step-tapered  undulator. 


Figure  12.  Positive  step-tapered  undulator  from  [3] . 


Figure  13.  Negative  step-tapered  undulator  from  [3] . 
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As  the  electron  beam' s  energy  changes  along  the 
undulator,  the  resonance  condition  is  altered.  By  tapering 
the  undulator,  we  introduce  the  following  phase  acceleration 
in  the  pendulum  equation  [7] : 

,  ,  AK 

S  =  -47tN - ^ (III.  2  0) 

\  +  K^  K 

where  AK/K  is  the  fractional  change  on  the  undulator 
parameter,  and  where  K  is  proportional  to  the  magnetic 
field . 

The  pendulum  equation  for  the  tapered  undulator  becomes 

[7]  : 


=  v  =  S  +  a  cos[i^ +  ^)  . 


(Ill  .21) 


The  optimal  extraction  from  a  normal  undulator  is  [7] : 


1 

2N 


but  in  a  tapered  undulator,  the  phase  acceleration 
an  important  role  in  the  extraction.  The  extraction 
case  is  [ 7 ]  : 


(III .22) 

5  plays 
for  this 


1^1 

n.  «  -L-L  .  ( 1 1 1 . 2  3 ) 

^  SttN 

As  the  phase  acceleration  (positive  or  negative) 
approaches  the  value  of  the  optical  field  amplitude,  the 
separatrix  becomes  smaller.  When  this  acceleration  is 
greater,  the  separatrix  vanishes,  and  bunching  does  not 
occur  because  the  electrons  are  not  trapped. 

It  is  obvious  that  at  bigger  phase  acceleration,  the 
extraction  will  be  bigger;  but  there  is  an  upper  limit.  This 
limit  is  set  by  the  fact  that  the  pendulum  equation  has  no 
solution  when  the  taper  phase  acceleration  is  bigger  than 
the  optical  field  amplitude;  however,  a  lower  limit  is 
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determined  by  the  maximum  deceleration  a  tapered  electron 
can  undergo  in  the  no  tapered  case.  Then,  we  can  state  the 
taper  range  as : 

4Va  <  <5’  <  |a| 

Figures  14  and  15  show  the  phase  space  evolution  of  a 
positive  and  negative  linearly  tapered  undulator.  In  the 
positive  taper  we  can  see  that  some  of  the  electrons  are 
trapped  inside  the  separatrix  and  they  bunch,  improving  the 
energy  exchange.  In  the  negative  taper  case,  we  can  see  that 
electrons  go  around  the  separatrix,  but  we  still  have 
bunching  and  good  extraction.  For  this  last  case,  the  taper 
will  work  well  only  when  electron  beam  starts  above 
resonance.  Optimal  taper  occurs  when  approximately  half  of 
the  electrons  are  trapped. 


Figure  14 . 


Phase  space  diagram  of  a  positive  linearly 
tapered  undulator  from  [3] . 
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Figure  15.  Phase  space  diagram  of  a  negative  linearly 

tapered  undulator  from  [3] . 

Figures  16  and  17  show  the  phase  space  evolution  of  a 
positive  and  negative  step-tapered  undulator.  There  we  can 
see  that  positive  taper  traps  electrons  like  the  positive 
linearly  tapered  case;  and  negative  tapered  electrons  are 
similar  to  the  negative  linearly  tapered  case. 


Figure  16. 


Phase  space  diagram  of  a  positive  step- 
tapered  undulator  from  [3] . 
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Figure  17 . 


Phase  space  diagram  of  a  negative  step- 
tapered  undulator  from  [3] . 


H.  DIFFRACTION 


Propagation  of  a  laser  beam  can  be  described  in  free 
space  by  Maxwell's  wave  equation  [7] 


C 


2  a,2 

C  Ot 


V 


a{x,t)  =  0  , 


(III  .24) 


where  a(x,t)e“^  = 


a(x,t) 


i(pa{x ,t)  ia 


e  is  the  optical  electric  field. 


For  laser  light,  the  beam  is  coherent  and  the  electric 
field  is  slow-varying  in  the  direction  of  propagation  over 
an  optical  wavelength  and  slow-varying  in  time  over  the 
optical  frequency.  In  the  last  equation,  a=  kz-cot  represents 
the  phase  of  the  carrier  wave,  and  the  optical  field  is 
represented  by  a,  as  defined  earlier. 

If  the  optical  field  is  constant,  the  electrical  field 
describes  a  simple  plane  wave  propagating  in  the  z 
direction,  containing  only  one  frequency  co .  If  the  wave's 
amplitude  or  phase  contain  spatial  and  time  dependence  the 
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laser  beam  will  have  a  finite-with  spectrum;  if  both 
dependences  are  slow  then  the  laser  will  have  a  narrow 
spectrum.  In  an  FEL,  the  amplitude  and  phase  are  slowly- 
varying,  so  we  can  make  the  following  assumptions: 
d\a\  ,,86  ,  d\a\  ,  ,  86 

—^«,o)\a\  ,  —<s,o)(b  ,  — ,  —<^k(i)  (III.  2  5) 

8t  '  '  5^  8z  '  '  az 

After  substitutions  we  get 


f  8^ a  - aa  2  ']  -2 

Via+  — :r  +  2ik - ka  -c 

8z^  8z 


8^ a  8a  ?  ^ 

— ^-2ico - ora  e  =0 

8t^  8t 


(III .26) 


Using  assumptions  (III. 25),  we  can  neglect  the  terms 

with  two  derivatives.  Multiplying  by  and  using  co=kc  we 

arrive  at 

(8  1  a  ^  ^ 

Vi+2zA: 1 - a(x,t)  =  Q  .  (III.  27) 

vaz  c  8t  Jj 

This  equation  can  be  simplified  further  by  introducing 
a  coordinate  transformation  u=z-ct  and  a  dimensionless  time 
T=ct/L  where  L  is  the  length  of  the  undulator.  Then,  the 
wave  equation  can  be  written  with  one  less  partial 
derivative , 

a(x,t)  =  0,  or  (III.  28) 


2ik  8 

Vi  + - 

L  8t 


a{x,t)  =  0  . 


(Ill  .29) 


This  equation  is  known  as  the  parabolic  wave  equation. 

The  first  term  in  this  equation  is  the  diffraction  term,  and 

its  magnitude  depends  on  the  transverse  area  of  the  laser 

beam,  the  wavelength  (X=2n/k) ,  and  the  undulator  length  L. 

The  Rayleigh  length  Zg  (distance  from  the  mode  waist  to 

where  the  area  of  the  optical  mode  doubles  due  to 

diffraction)  relates  the  optical  mode  waist  radius  to 
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the  wavelength  of  the  light  by  ZQ=7tWQlX.  In  dimensionless 
units  Zq=Z^IL  .  Based  on  diffraction  along  the  undulator,  so 
Z^'^L,  the  characteristic  mode  waist  radius  is 
4UJ  7t  =-sj2LI k  .  If  the  optical  mode  radius  is  much  larger 

than  \J LA/ 71  ,  then  the  diffraction  has  a  small  effect  over 
the  undulator.  If  it  is  smaller,  diffraction  has  a  big 
effect  on  propagation,  changing  both  the  amplitude  and  phase 
of  the  optical  field. 

It  is  important  to  mention  that  this  description  of 
diffraction  is  calculated  for  laser  light  in  free  space.  For 
an  FEL  in  the  amplifier  configuration,  diffraction  is 
altered  by  the  electron  interaction.  Inside  the  undulator, 
the  wave  equation  is  inhomogeneous,  and  the  Gaussian 
distribution  of  the  optical  field  changes;  the  laser  beam  is 
affected  by  the  electrons  that  try  to  focus  the  laser  beam. 
Then  we  have  to  consider  the  last  term  of  the  FEL  wave 
equation  <-je“"  >  .  This  new  term  affects  the  diffraction  of 
the  optical  field. 

I.  SHORT  PULSES  AND  SLIPPAGE 

Eor  energy  exchange  between  the  electrons  and  optical 
field,  they  must  overlap  inside  the  undulator.  During  their 
pass  through  the  wiggler,  electrons  and  photons  have 
different  velocities.  Photons  travel  at  speed  of  light,  and 
relativistic  electrons  move  close  to  it  but  slower.  The  term 
used  to  describe  the  overlap  of  the  optical  and  electron 
pulses  is  known  as  "slippage";  according  to  the  resonance 
condition,  the  electrons  drop  behind  the  light  by  a  slippage 
distance  NX  over  the  length  of  the  undulator.  Eigure  18 
shows  the  slippage  concept.  Red  pulse  is  the  electron  pulse 
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and  blue  pulse  is  the  optical  pulse;  the  left  overlap 
represents  the  position  at  the  beginning  of  the  undulator, 
and  the  right  overlap  is  at  the  end  of  the  undulator. 


Figure  18.  Electron  and  light  slippage  from  [3]  . 


We  have  short  pulses  when  the  micropulse  length  ^  is 

on  the  order  of  or  less  than  the  slippage  distance  NX.  In 
order  to  get  the  required  energy  transfer,  the  short 
electron  pulses  and  optical  pulses  must  be  synchronized. 
That  process  develops  in  stages,  and  electron  bunching  must 
take  place  in  the  right  moment.  At  the  beginning  of  the 
undulator  the  electrons  and  light  overlap,  but  there  is  no 
bunching  and  therefore  no  gain.  At  the  middle,  electrons 
during  their  bunching  process  slightly  amplify  the  light 
pulse.  At  the  end,  bunched  electrons  amplify  the  trailing 
part  of  the  optical  pulse.  As  a  result,  the  optical  pulse 


29 


centroid  travels  slower  than  speed  of  light;  this  process  is 
called  "optical  lethargy".  Figure  19  shows  this  process. 


undulator 


Figure  19.  Optical  lethargy  from  [3] 
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IV.  SIMULATIONS  OF  FREE  ELECTRON  LASERS  IN 
AMPLIFIER  CONFIGURATION 


In  this  chapter,  several  free  electron  laser  amplifiers 
were  explored  with  computer  simulations,  using  the  NFS  Apple 
Xserve  cluster  with  64  nodes/128  processors.  On  this  device, 
a  simulation  program  was  run  that  solved  the  coupled 
pendulum  and  FEL  wave  equations  in  order  to  understand  the 
four  dimensional  behavior  of  the  high  power  FEL  amplifier. 
The  simulation  program  required  dimensionless  input 
parameters,  which  make  the  results  applicable  to  a  broad 
range  of  EELs . 

A.  DESCRIPTION  OF  THE  INITIAL  PARAMETERS 

The  first  step  in  this  process  is  to  establish  the 
initial  parameters  that  describe  the  system  to  be  evaluated. 
Using  these  parameters,  the  dimensionless  input  variables 
are  calculated  in  order  to  run  the  simulation. 

We  can  classify  these  parameters  as  primary  and 
secondary.  Secondary  parameters  are  coupled  to  the  primary 
parameters  and  they  have  to  be  calculated  using  several 
formulas  that  will  be  described  in  this  work. 

The  desired  characteristics  of  the  device  were 
introduced  into  an  Excel  spread  sheet,  which  calculated  the 
secondary  parameters  using  formulas,  and  the  dimensionless 
input  parameters  for  the  simulation  program.  Tables  1,  2  and 
3  show  the  description  of  the  dimensional  parameters  that 
were  used  and  the  formulae  used  to  obtain  the  secondary 
dimensional  parameters. 
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H undulator  parameters 

/Ig 

Undulator  period 

N 

Number  of  periods 

g 

Undulator  gap 

B 

Magnetic  field  (rms) 

K 

Undulator  parameter 

—  SBXg/ 

/  2  Tzmc^ 

L 

Undulator  length 

II 

Table  1 .  Undulator  parameters  used  to  describe  an  FEL  in  a 

simulation . 


ELECTRON  BEAM  PARAMETERS 

E, 

Beam  kinetic  energy 

(It 

Bunch  charge 

h 

Beam  radius  (rms) 

h 

Pulse  duration,  EWHM 

Q 

Pulse  repetition 
frequency 

4 

Pulse  length,  EWHM 

4  =  he 

r 

Lorentz  factor 

^  jE,+t.5nMeV) 

/ 

OAUMeV 

^  Peak 

Peak  current 

I  = 

Peak  /  f 

/ 

Pvg 

Average  current 

II 

ID 

Normalized  rms  emittance  \ 

1 

Longitudinal  emittance 

A/ 

r 

Beam  energy  spread  (rms) 

^  _  2.35^;/ 
y  /^Eih){t^) 

AO 

Beam  angular  spread 
(rms) 

A6»  =  ^«/ 

/Yh 

Pt 

Beam  average  power 

P  -  F  f 

^  Avg 

Table  2 .  Electron  beam  parameters  used  to  describe  an  EEL 

in  a  simulation. 
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■  OPTICAL  PARAMETERS  ' 

P,n 

Seed  laser  power 

Seed  pulse  duration 

s 

Distance  to  First  Optic 

Rayleigh  length 

Z 

^0  “  /Z 

Optical  wavelength 

/ly^ 

w. 

Mode  waist  radius,  1/e 

II 

Table  3.  Optical  parameters  used  to  describe  an  FEL  in  a 

simulation . 


B.  SIMULATIONS  OF  THE  FIRST  FEL 

This  first  laser  was  intended  to  have  a  power  of  100  kW 
with  a  wavelength  of  1.6  pm  in  an  amplifier  configuration. 
The  undulator  has  120  periods,  and  its  length  is  324  cm. 
Tables  4  through  7  show  these  parameters. 


I  ELECTRON  BEAM  PARAMETERS 

E, 

Beam  energy  (MeV) 

80 

9b 

Bunch  charge  (nC) 

0.2 

h 

Beam  radius,  (mm) 

0.13 

h 

Pulse  duration,  FWHM  (ps) 

1.0 

Q 

Pulse  rep  frequency  (MHz) 

703 

4 

Pulse  length,  FWHM  (cm) 

0.03 

y 

Lorentz  factor 

158 

^ Peak 

Peak  current  (A) 

200 

^  Avg 

Average  current  (mA) 

141 

Normalized  rms  emittance  (mm  mrad) 

3.5 

^1 

Longitudinal  emittance  (keV  ps) 

36 

<] 

Beam  energy  spread  (%) 

0.11 

A0 

Beam  angular  spread  (mrad) 

0.17 

Pb 

Beam  average  power  (MW) 

11 

Initial  electron  beam  parameters  of  the  first 
simulated  FEL  device. 

33 


Table  4 . 


UNDULATOR  PARAMETERS 


Ag  Undulator  period  (cm)  2.7 
N  Number  of  periods  120 
g  Undulator  gap  (cm)  1.00 
K  Undulator  parameter,  rms  1.39 
L  Undulator  length  (cm)  324 


Table  5.  Initial  undulator  parameters  of  the  first 

simulated  FEL  device. 


Table  6.  Initial  optical  parameters  of  the  first  simulated 

FEL  device. 


Table  7.  Dimensionless  parameters  of  the  first  simulated 

FEL  device. 
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The  first  parameter  to  analyze  was  the  taper  start  time 
.  This  time  is  established  as  a  dimensionless  parameter 
using  a  scale  from  0  to  1,  with  r^=0  corresponding  to  the 
beginning  of  the  undulator  and  ^^=1  the  end  of  he 
undulator.  The  taper  start  time  was  simulated  for  several 
values  from  0.4  to  0.8  with  intervals  of  0.1;  using  a  taper 
rate  of  5=40n.  The  start  time  with  the  greatest  extraction 
was  at  r^=0.7  with  an  extraction  g  =  0.92%.  Figure  20  shows 
a  graph  of  the  obtained  extractions  versus  the  simulated 
taper  start  times  for  5=40n. 


Figure  20.  Graph  of  extraction  vs.  taper  start  times. 

The  full  simulation  output  for  the  most  efficient  taper 
start  time  ( =  0.7 )  at  5=40n  is  shown  on  Figure  21.  There  we 
can  see  the  power  P(t)  growing  to  saturation  at  the  end  of 
the  undulator  (t=1),  with  about  half  the  electrons  trapped 
in  phase  space  (^,v)  as  shown  on  phase  space  diagram.  Also 
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we  can  see  a  typical  diffraction  pattern  of  the  beam  in  x 
and  y.  Slippage  produces  a  narrow  optical  pulse  in  the  z 
direction . 


Simulation  run  time:  39.5  sec  Fri  May  4  11:31:22  2007 

nx=200,  nz=33,  nt=2  5G,  np=100000,  Wp=12,  seed=7,  v4jins=33,  ebins=16 


2  -2 


2  -4 


4-13 


13 


0.0  X  1.0  0m7-2x  2-5x  5 

F=0.44,  G=1.14e+03,  ti=0.00921,  ^7/7=0. 0408 
THEORY:  G(jF»l)=5.  14e+04 
THEORY:  r| ( jF»l)=0.  0125 

a  =141  at  X  =0.70,  THEORY:  a  =164  at  x  =0.57 
s  s  s  s 

M^=0.047,  c^(0,0)=0.25,  (0, 2) =0. 11,  c^(2,0)=0.11 


Figure  2 1 . 


Full  simulation  output  for  taper  start  time 

0.7. 
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Once  the  most  efficient  value  of  the  taper  start  time 
(r^=0.7)  was  found,  it  was  used  to  run  a  set  of  simulations 
to  analyze  the  dimensionless  linear  taper  rate  (5)  . 
Simulations  were  run  for  5  values  from  0  to  80n  at  intervals 
of  lOn.  The  most  efficient  rate  was  at  5  =  40n  with  an 

extraction  r]  =  0.92%.  Figure  22  shows  a  graph  of  the 

obtained  extractions  versus  simulated  taper  rates. 


Figure  22.  Graph  of  extraction  vs.  linear  taper  rates. 


The  full  simulation  output  for  the  most  efficient  taper 
start  time  ( =  0.7 )  is  shown  in  Figure  21  ( coincidently  the 
taper  rate  used  in  that  simulation  was  the  most  efficient) . 
The  results  are  the  same  as  the  previous  simulation.  Since 
the  extraction  is  not  enough  to  achieve  100  kW,  it  is 
necessary  to  change  other  parameters  in  order  to  reach  the 
goal . 


37 


In  order  to  look  for  a  better  extraction,  a  longer 
undulator  was  tried,  increasing  the  number  of  periods  to 
140.  Since  the  length  of  the  undulator  changed  to  378  cm,  it 
was  necessary  to  calculate  the  dimensionless  input 
parameters  again.  Table  8  shows  only  the  parameters  that 
changed  for  this  case. 


I  ELECTRON  BEAM  PARAMETERS 

^  < 

Beam 

Beam 

radius , 
angular 

(mm) 

spread  (mrad) 

0.14 

0.15 

Table  8.  New  parameters  for  N=140. 


Using  a  taper  start  time  r^=0.5,  a  set  of  simulations 
explored  the  effect  of  the  dimensionless  taper  rate  (5) . 
Simulations  were  run  for  5  values  from  0  to  80n  at  intervals 
of  lOn.  The  most  efficient  rate  was  at  5  =  30n  with  an 
extraction  g  =  1.13%.  Figure  23  shows  a  graph  of  obtained 
extractions  versus  the  simulated  taper  rates. 
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Dimensionless  Linear  Taper  Rate  in  Multiples  of  pi 

Figure  23.  Graph  of  extraction  vs.  linear  taper  rates 

for  N=140. 


With  this  longer  undulator  we  can  see  that  the 
extraction  is  bigger  even  with  a  smaller  taper  rate. 
Diffraction  patterns  are  very  similar  to  the  previous  case 
and  a  higher  power  reaches  saturation  at  the  end  of  the 
undulator.  On  Figure  24  we  can  see  that  electron  bunching  is 
quite  similar  to  the  previous  case,  but  the  optical  field 
evolved  into  a  combination  of  higher-order  modes. 


-11/2  ^  3ii/2  -2x2 


Figure  24.  Electron  bunching  and  optical  field  for 

linear  taper  rate  5=30n. 
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Continuing  with  this  exploration,  and  using  the  optimal 
taper  rate  found  for  this  FEL  (5=30n)  and  a  taper  start  time 
r^=0.6,  this  device  was  analyzed  by  changing  the  electron 
beam  pulse  duration  ( )  in  a  range  from  0.4  to  1.6  ps  in 
0.2  ps  steps.  It  is  important  to  mention  that  any  change  in 
this  parameter  changes  the  input  parameters  of  the 
simulation,  so  a  recalculation  of  dimensionless  parameters 
was  made  for  every  value  of  .  The  largest  extraction  (g  = 
1.18%)  was  found  at  ^^=0.8ps.  Figure  25  shows  a  graph  of  the 
obtained  extractions  versus  the  simulated  electron  beam 
pulse  duration. 


Figure  25.  Graph  of  extraction  vs.  electron  beam  pulse 

duration . 
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The  results  at  the  best  pulse  duration  (^j=0.8ps)  are 
similar  to  the  last  case,  but  this  time  it  did  not  have  many 
higher-order  optical  modes. 

The  next  step  is  to  change  other  parameters  to  obtain  a 
better  extraction.  This  time  the  bunch  charge  was  changed 
from  0.2  to  0.4  nC,  with  the  pulse  repetition  frequency  cut 
in  half  (Q  =  350MHz)  to  keep  the  average  current  the  same 
(140  mA)  .  This  change  affects  some  of  the  input  parameters 
to  the  simulation  program,  so,  it  is  necessary  to 
recalculate  the  dimensionless  parameters  to  create  another 
input  file.  Table  9  shows  the  initial  parameters  that 
changed  for  this  new  case  with  respect  to  the  last  case  upon 
changing  the  bunch  charge  value. 


Table  9.  Parameters  that  changed  for  micropulse  charge  of 

0.4 
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Using  these  new  parameters  the  input  file  was  updated 
for  the  simulation,  and  several  runs  tried  to  obtain  optimal 
outputs  for  the  same  ranges  of  taper  values. 

Again  the  first  parameter  to  analyze  was  the  taper 
start  time.  The  taper  start  time  was  simulated  for  several 
values  from  0.4  to  0.8  at  intervals  of  0.1;  using  a  linear 
taper  rate  of  40n.  The  optimum  start  time  was  at  0.5  with  an 
extraction  g  =  2.06%.  Figure  27  shows  a  graph  of  the 
obtained  extractions  versus  the  simulated  taper  times. 


Figure  26.  Extraction  vs.  dimensionless  taper  start  time 

for  micropulse  charge  of  0.4  nC . 

The  results  for  the  optimum  taper  time  (0.5)  are  very 
similar  to  the  last  case,  but  this  there  was  slightly  better 
electron  bunching,  resulting  in  a  better  extraction.  As  a 
result  and  a  higher  power  was  delivered  when  saturation  was 
reached  at  the  end  of  the  undulator. 
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Again,  once  the  optimum  value  of  the  taper  start  time 
(0.5)  was  found,  it  was  used  to  run  a  set  of  simulations  to 
analyze  the  dimensionless  taper  rate  (5)  .  Simulations  were 
run  for  6  values  from  0  to  80n  at  intervals  of  lOn.  The 
optimum  rate  was  at  5  =  60n  and  5  =  70n  with  an  extraction  p 
=  2.32%.  Figure  29  shows  a  graph  of  the  obtained  extractions 
versus  the  simulated  taper  rates. 


Figure  27.  Extraction  vs.  dimensionless  taper  rates  for 

micropulse  charge  of  0.4  nC . 

The  results  for  the  taper  rate  5  =  60n  are  shown  in 
Figure  28.  Extraction  is  much  better  in  this  case,  but  the 
taper  rate  must  be  large  in  order  to  reach  this  level  of 
extraction.  As  in  the  last  case,  a  higher  power  was 
delivered  when  saturation  was  reached  at  the  end  of  the 
undulator.  It  may  be  possible  to  reach  these  levels  changing 
other  parameters . 
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simulation  run  time:  40.7  sec  Tue  Hay  8  13:41:02  2007 

rDC=200,  nz=33,  nt=2  56,  r]p=100000,  Wp=12,  seed=7,  vftiins=33,  ebins=16 


*****  pel  amplifier 

4D  SIMULATION  *****  L 

|a| 

max 

j=9137  aj^=0.1 

ay=0.1 

a^=0.42 

a^=0.42  a^=1.3 
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II 

o 
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5  =0ii  X  =0 . 6 

a  =9.7 

Q 

II 

t  =0  z  =0.09 

X  =3 

0  =0 

s  s 

o 

a 

w  o 

m 

m 

1 

X 

0 


-2  X  2-2  Y  2-2  z  2-25  v  25 


0.0  X  1.0  Om7-2x  2-5x  5 

F=0.44,  G=2.7e+03,  ti=0.0232,  A7/7=0.078  5 
THEORY:  G(jr»l)=l.  23e+04 
THEORY:  ti (jF»l)=0.  0214 

a  =329  at  x  =0.60,  THEORY:  a  =321  at  x  =0.42 
s  s  s  s 

M^=0.025,  c^(0,0)=0.24,  (0, 2) =0. 11,  c^(2,0)=0.11 


Figure  28.  Results  for  micropulse  charge  of  0.4  nC  and 

taper  rate  5=60n. 
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The  last  explored  parameter  was  the  emittance;  this 
parameter  is  mainly  determined  by  the  injector.  Emittance 
changes  have  a  direct  effect  on  the  initial  parameters;  and 
any  change  of  this  parameter  needs  to  be  applied  to  the 
spread  sheet  in  order  to  calculate  the  new  input  parameters. 

For  this  simulation,  emittance  was  adjusted  to  5.1  mm 
mrad,  and  Table  10  shows  the  initial  parameters  that  changed 
for  this  new  case  with  respect  to  the  last  case  upon 
changing  emittance. 


Normalized  rms  emittance  (mm  mrad)  5.1 

A0  Beam  angular  spread  (mrad)  0.23 


1  OPTICAL  PARAMETERS 

Wg  Mode  waist  radius,  1/e  (mm) 

0.40 

Table  10.  Parameters  that  changed  after  changing  the 

emittance . 


With  these  parameters  and  using  a  taper  start  time 
r^=0.5,  it  was  run  a  set  of  simulations  to  explore  the 
dimensionless  taper  rate  (5)  .  Simulations  were  run  for  5 
values  from  30n  to  80n  at  intervals  of  lOn.  The  optimum 
taper  rate  was  at  5  =  60n  with  an  extraction  g  =  2.15%. 
Figure  29  shows  a  graph  of  obtained  extractions  versus  the 
simulated  taper  rates. 
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2.5 


Dimensionless  Linear  Taper  Rate  in  Multiples  of  pi 


Figure  29.  Extraction  vs.  dimensionless  taper  rates  for 

emittance=5 . 1mm  mrad. 

This  graph  shows  that  with  this  value  of  emittance, 
taper  rate  variation  over  this  range  does  not  have  a  big 

effect  on  the  extraction.  For  5  =  60n,  the  results  of  the 
simulation  are  very  similar  to  the  last  case.  As  in  the 

last  case,  a  higher  power  was  delivered  when  saturation  was 
reached  at  the  end  of  the  undulator. 

C.  SIMULATIONS  OF  THE  SECOND  FEL 

This  second  set  of  simulations  explores  an  FEL 
amplifier  with  a  power  of  1  megawatt  and  an  optical 

wavelength  of  1.6  microns.  The  same  procedure  is  used  for 

this  laser,  but  this  time  the  number  of  periods  of  the 
undulator  was  analyzed  for  several  values  in  order  to  find 
the  optimum.  The  general  input  parameters  of  this  device  are 
shown  on  Tables  11,  12,  13  and  14. 
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ELECTRON  BEAM  PARAMETERS 


E, 

Beam  energy  (MeV) 

75.5 

(It 

Bunch  charge  (nC) 

1.5 

Beam  radius,  x  (mm) 

0.31 

hy 

Beam  radius,  y  (mm) 

0.18 

h 

Pulse  duration,  FWHM  (ps) 

1 . 1 

a 

Pulse  rep  frequency  (MHz) 

703 

h 

Pulse  length,  FWHM  (cm) 

0.034 

r 

Lorentz  factor 

149 

^  Peak 

Peak  current  (A) 

1333 

^ Avg 

Average  current  (mA) 

1055 

^nx 

Normalized  rms  emittance,  x  (mm 

mrad)  15.0 

^ny 

Normalized  rms  emittance,  y  (mm 

mrad)  5 . 0 

Longitudinal  emittance  (keV  ps) 

70 

Af 

7 

Beam  energy  spread  (%) 

0.19 

Beam  angular  spread,  x  (mrad) 

0.32 

Beam  angular  spread,  y  (mrad) 

0.18 

Pt 

Beam  average  power  (MW) 

80 

Table  11.  Initial  electron  beam  parameters  of  the  second 

simulated  FEL. 


I  UNDULATOR  PARAMETERS 

Aq 

Undulator  period  (cm) 

2 . 62 

N 

Number  of  periods 

150 

g 

Undulator  gap  (cm) 

1.00 

K 

Undulator  parameter,  rms 

1.30 

L 

Undulator  length  (cm) 

393 

Table  12.  Initial  undulator  parameters  of  the  second 

simulated  FEL  device. 
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OPTICAL  PARAMETERS 


Pi„  Seed  laser  power  (W)  100 

T.^  Seed  pulse  duration,  FWHM  (ps)  2.3 

S  Distance  to  First  Optic  (cm)  1000 

Zg  Rayleigh  length  (cm)  105 

X  Optical  wavelength  (microns)  1.59 

Wq  Mode  waist  radius,  1/e  (mm)  0.73 

Table  13.  Initial  optical  parameters  of  the  second  simulated 

FEL  device. 


DIMENSIONLESS  PARAMETERS 


j  Normalized  current  density,  linear  undulator  20288 
cr^  Normalized  beam  radius,  x  0.22 

cr^.  Normalized  beam  radius,  y  0.13 

(J,^  Normalized  beam  angular  spread,  x  0.8  9 

(Jty  Normalized  beam  angular  spread,  y  0.52 

<J^g  Phase  velocity  spread  due  to  energy  spread  3.7 

cr^  Normalized  pulse  length  1.4 

flg  Normalized  initial  optical  field  amplitude  4.3 

cr^  Normalized  initial  optical  pulse  length  2.8 

Zq  Normalized  Rayleigh  length  0.27 

Wg  Normalized  mode  waist  radius  0.52 

co^  Betatron  oscillation  frequency  8.2 


Table  14.  Dimensionless  parameters  of  the  second  simulated 

FEL  device. 
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For  the  first  simulation,  it  was  necessary  to  update 
the  dimensionless  input  file  for  each  value  of  N.  The  number 
of  periods  was  analyzed  from  N=100  through  N=200  at 
intervals  of  10,  obtaining  the  highest  extraction  g  =  2.21% 
at  N=140.  Figure  30  shows  a  graph  of  the  extraction  versus 
the  number  of  periods  of  the  undulator. 


Figure  30.  Extraction  vs.  number  of  periods  of  the 

second  FEL. 

Figure  31  shows  the  results  of  the  simulation  for  this 
set  of  parameters  with  N=140,  which  is  the  optimum  value. we 
can  see  the  power  P(t)  growing  to  saturation  at  the  end  of 
the  undulator  (t=1)  with  an  intermediate  bump  corresponding 
to  a  synchrotron  oscillation,  and  about  half  the  electrons 
trapped  in  phase  space  (^,v)  as  shown  on  phase  space 
diagram. 


49 


Simulation  run  time:  39.9  sec  Tue  Aug  7  13:05:59  2007 

nx=200,  nz=33,  nt=2  56,  1^5=100000,  Wp=12,  seed=7,  ’wbins=33,  ebins=20 


139 


-139 


P(x)  2.87xl0^|-’t/2  ^  3ji/2  -3 


0.0 


1.0 


-3 


F=0.42,  G=2.03e+04,  ti=0.0221,  A7/ic0.09  54 
THEOror:  G(jF»l)=2.  65e+04 
THEOrar:  Ti(jF»l)=0.  0179 

a  =535  at  x  =0.50,  THEORY:  a  =371  at  x  =0.48 
s  s  s  s 

M^=0.058,  c^(0,0)=0.33,  (2, 0) =0. 13,  c^(0,2)=0.12 


Figure  31. 


Results  for  N=140  of  the  second  FEL. 
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Using  the  undulator  with  N=140,  the  taper  start  time 
was  simulated  for  several  values  from  0.0  to  0.9  at 
intervals  of  0.1  with  a  taper  rate  of  40n.  The  optimum  start 
time  (greatest  extraction)  was  at  0.5  with  an  extraction  g  = 
2.21%.  Figure  32  shows  a  graph  of  the  obtained  extractions 
versus  the  simulated  taper  times. 


Figure  32.  Extraction  vs.  dimensionless  taper  start  time 

of  the  second  FEL. 

The  results  for  the  best  taper  time  (0.5)  are  shown  in 
Figure  31. 

The  taper  start  time  of  0.5  was  used  to  run  a  set  of 
simulations  to  analyze  the  dimensionless  taper  rate  (5)  . 
Simulations  were  run  for  5  values  from  0  to  150n  at 
intervals  of  lOn.  The  optimum  rate  was  at  5  =  lOOn  with  an 
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Figure  33  shows  a  graph  of  the 


extraction  g  =  3.55%. 

obtained  extractions  versus  the  taper  rates. 


Figure  33.  Extraction  vs.  dimensionless  linear  taper 

rates  of  the  second  FEL. 

The  results  for  the  optimum  taper  rate  5  =  lOOn  are 
shown  in  Figure  34.  Extraction  is  much  better  in  this  case. 
As  in  the  last  case,  a  higher  power  was  delivered  when 
saturation  was  reached  at  the  end  of  the  undulator.  Less 
than  half  of  the  electrons  bunched  inside  of  the  separatrix 
as  shown  on  the  phase  space  diagram. 
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Simulation  run  time:  39.3  sec  Ned  Aug  8  12:23:38  2007 

nx=200,  nz=33,  nt=256,  i^=100000,  Wp=12,  seed=7,  vA>ins=33,  ebins=20 


*****  fel  amplifier  4d  simulation  *****  L - — - 1 

0  |a  I  max 
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Oy=0.l3 
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fl 

-3 


2  58 


-258 


3  -3 


3  -4 


P(x) 

4.52x10^ 

G(x) 

(;,v) 


. 


4-13  V  13 
a(x,y,0) 


-x/2  t,  37t/2  -3 


0.0  X  1.00  m  7  -3  x  3  -8  x  8 

F=0.42,  G=3.2e+04,  ti=0.0355,  A7/7=0.139 
THEORY:  G(jF»l)=6.  Ole+04 
THEORy:  r|  ( jF»l)=0.  0446 

a  =489  at  x  =0.50,  THEORY:  a  =371  at  x  =0.48 
s  s  s  s 

M^=0.094,  c^(0,0)=0.35,  c^(2,0)=0. 14,  c^(0,2)=0.12 


Figure  34.  Results  for  5  =  lOOn  of  the  second  FEL. 
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D. 


SIMULATIONS  OF  THE  THIRD  FEL 


This  third  set  of  simulations  explores  an  FEL  amplifier 
with  a  power  of  3  megawatts  and  an  optical  wavelength  of  1.6 
microns.  In  this  case,  more  realistic  electron  beam 
parameters  were  generated  by  an  external  computer  program 
(Parmela) .  The  number  of  periods  of  the  undulator  was 
adjusted  in  order  to  obtain  the  most  extraction.  The  general 
input  parameters  of  this  device  are  shown  on  Tables  15 
through  1 8 . 


ELECTRON  BEAM  PARAMETERS 


Beam  energy  (MeV)  75.5 

qi,  Bunch  charge  (nC)  1.5 

Beam  radius,  x  (mm)  0.2  8 

Beam  radius,  y  (mm)  0.18 

Pulse  duration,  FWHM  (ps)  0.88 

Q  Pulse  rep  frequency  (MHz)  703 

4  Pulse  length,  FWHM  (cm)  0.026 

/  Lorentz  factor  149 

Peak  current  (A)  17  05 

Average  current  (mA)  1055 

Normalized  rms  emittance,  x  (mm  mrad)  14.0 

s„y  Normalized  rms  emittance,  y  (mm  mrad)  5.6 

Longitudinal  emittance  (keV  ps)  60 

A/ 

-  Beam  energy  spread  (%)  0.21 

r 

Ad^  Beam  angular  spread,  x  (mrad)  0.33 

A0y  Beam  angular  spread,  y  (mrad)  0.21 

Pi,  Beam  average  power  (MW)  80 


Table  15.  Initial  electron  beam  parameters  of  the  third 

simulated  FEL. 
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UNDULATOR  PARAMETERS 


Aq 

Undulator 

period  (cm) 

2 . 62 

N 

Number  of 

periods 

130 

g 

Undulator 

gap  (cm) 

1.00 

K 

Undulator 

parameter,  rms 

1.30 

L 

Undulator 

length  (cm) 

341 

Table  16.  Initial  undulator  parameters  of  the  third 

simulated  FEL  device. 


Table  17.  Initial  optical  parameters  of  the  third  simulated 

FEL  device. 


Table  18.  Dimensionless  parameters  of  the  third  simulated 

FEL  device. 
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For  the  first  simulation,  it  was  necessary  to  update 
the  dimensionless  input  file  each  time  the  number  of  periods 
was  changed.  The  number  of  periods  was  varied  from  N=100 
through  N=200  in  intervals  of  10,  obtaining  the  highest 
extraction  of  ri=3.15%  at  N=130.  Figure  35  shows  a  graph  of 
the  extraction  versus  the  number  of  periods  of  the 
undulator.  We  also  considered  rotating  the  undulator  90°  by 
switching  the  x  and  y  components  of  the  normalized  beam 
radius  and  angular  spread.  This  change  did  not  make  any 
significant  difference  in  the  results. 


Figure  35.  Extraction  vs.  number  of  periods  of  the  third 

FEL. 
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Figure  38  shows  the  simulation  results  for  this  set  of 
parameters  with  N=130,  which  is  the  optimal  value.  We  can 
see  the  power  P(t)  growing  to  saturation  at  the  end  of  the 
undulator  (t=1)  with  an  intermediate  synchrotron 
oscillation,  and  about  half  of  the  electrons  bunched  in 
phase  space  (^,v)  as  shown  on  phase  space  diagram. 

Using  the  optimum  undulator  (N=130) ,  the  taper  start 
time  was  simulated  for  several  values  from  0.2  to  0.8  at 
intervals  of  0.1  using  a  taper  rate  of  80n.  The  optimal 
start  time  (greatest  extraction)  was  at  0.5  with  an 
extraction  g  =  3.14%.  Figure  36  shows  a  graph  of  the 
extractions  versus  the  taper  start  times. 


Figure  36.  Extraction  vs.  dimensionless  taper  start  time 

of  the  third  FEL. 
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The  results  for  the  optimal  taper  start  time  (0.5)  are 
shown  in  Figure  38. 

Again,  once  the  optimal  value  of  the  taper  start  time 
(0.5)  was  found,  it  was  used  to  run  a  set  of  simulations 
varying  the  dimensionless  linear  taper  rate  (5) .  Simulations 
were  run  for  5  values  from  20n  to  llOn  at  intervals  of  lOn. 
The  most  efficient  rates  were  at  5  =  80n  and  5  =  90n  with  an 
extraction  g  =  3.14%.  Figure  37  shows  a  graph  of  the 
obtained  extractions  versus  the  taper  rates.  Beyond  the  peak 
it  is  roughly  flat. 


Figure  37.  Extraction  vs.  dimensionless  taper  rates  of 

the  third  FEL. 

The  results  for  the  taper  rate  5  =  80n  are  shown  in 
Figure  38. 
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simulation  run  time:  31.9  sec  Tue  Aug  21  10:09:38  2007 

nx=300,  nz=33,  nt=2  5G,  i5)=30000,  ^3=16,  seed=7,  vAiins=38,  eblns=15 
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Figure  38.  Results  for  N=130  of  the  third  FEL. 
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V.  CONCLUSION 


In  this  thesis,  we  have  discussed  the  components  and 
basic  theory  of  a  Free  Electron  Laser  (EEL)  in  order  to 
describe  its  behavior. 

We  have  explored,  by  means  of  simulations,  three 
designs  of  Eree  Electron  Lasers  in  the  amplifier 
configuration,  using  different  parameters  to  establish  a 
desired  power  output  and  optical  wavelength.  It  has  been 
shown  that  undulator  taper  has  an  important  effect 
increasing  the  single  pass  extraction. 

Eor  the  low  power  (~100kW)  EEL  examples,  extraction 
between  1%  and  2%  was  found.  Even  after  optimizing  the  taper 
parameters,  the  extraction  did  not  improve  much.  We  tried  to 
obtain  higher  extraction  by  changing  electron  beam 
parameters,  such  as  pulse  duration,  bunch  charge,  and 
emittance,  but  found  no  significant  improvement.  We  conclude 
that  these  three  parameters  did  not  have  a  large  effect  on 
the  final  extraction  of  the  system. 

In  higher  power  (~3MW)  EEL  examples,  it  was  found  to  be 
easier  to  obtain  higher  extraction.  Eor  high  power  EEL 
examples,  we  have  obtained  extractions  above  3%,  compared  to 
1%  or  2%  in  low  power  devices. 

Euture  investigation  should  focus  on  taper  parameters 
and  undulator  characteristics,  as  an  important  way  to 
improve  EEL  designs. 
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